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The existence of drag reduction by polymer additives, well established for wall-bounded turbulent 
flows, is controversial in homogeneous, isotropic turbulence. To settle this controversy we carry out 
a high-resolution direct numerical simulation (DNS) of decaying, homogeneous, isotropic turbulence 
with polymer additives. Our study reveals clear manifestations of drag-reduction-type phenomena: 
On the addition of polymers to the turbulent fluid we obtain a reduction in the energy dissipation 
rate, a significant modiflcation of the fluid energy spectrum especially in the deep-dissipation range, 
a suppression of small-scale intermittency, and a decrease in small-scale vorticity filaments. 
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The dramatic reduction of drag by the addition of 
small concentrations of polymers to a turbulent fluid con- 
tinues to engage the attention of engineers and physicists. 
Significant advances have been made in understanding 
drag reduction both experimentally 0, and theo- 
retically 0, 0, 0, i'^ channel flows or the Kolmogorov 
flow [Sj. However, the existence of drag-reduction-typc 
phenomena in turbulent flows that are homogeneous and 
isotropic [1 [11 [II [11 [11 [11 [11 fll remains 
controversial. Some experimental |15ij l(l|l7l Il8l | , numer- 
ical 0, 0, 0, Q| , and theoretical |9lll0l| studies have 
suggested that drag reduction should occur even in ho- 
mogeneous, isotropic turbulence; but other studies have 
refuted this claim [l9| . 

To settle this controversy we have initiated an exten- 
sive direct numerical simulation (DNS) of decaying, ho- 
mogeneous, isotropic turbulence in the presence of poly- 
mer additives. We monitor the decay of turbulence from 
initial states in which the kinetic energy of the fluid is 
concentrated at small wave vectors; this energy then cas- 
cades down to large wave vectors where it is dissipated 
by viscous effects; the energy-dissipation rate e attains 
a maximum at t„i, roughly the time at which the cas- 
cade is completed. A recent shell-model study ^1 has 
suggested that this peak in e can be used to quantify 
drag reduction by polymer additives. Since shell mod- 
els are far too simple to capture the complexities of real 
flows, we have studied decaying turbulence in the Navier- 
Stokes (NS) equation coupled to the Finitely Extensible 
Nonlinear Elastic Pcterhn (FENE-P) model for poly- 
mers. Our study, designed specifically to uncover drag- 
reduction-typc phenomena, shows that the position of 
the maximum in e depends only mildly on the polymer 
concentration c; however, the value of e at this maxi- 
mum falls as c increases. We use this decrease of e to de- 
fine the percentage drag (or dissipation) reduction DR in 
decaying homogeneous, isotropic turbulence; we also ex- 
plore other accompanying physical effects and show that 
they are in qualitative accord with drag-reduction exper- 



iments |la|2l|: In particular, DR increases with c (upto 
25% in one of our simulations). For small values of c 
the energy spectrum of the fluid is modified appreciably 
only in the dissipation range; however, this suffices to 
yield significant drag reduction. We show that vorticity 
filaments and intermittency are reduced at small spatial 
scales and that the extension of the polymers decreases 
as c increases. 

The NS and FENE-P (henceforth NSP) equations are 

Dtu = iyV^u+-i^V.[f{rp)C]-Vp; (1) 

TP 

DtC ^ C.(Vu) + (Vu)^.C- (2) 

Here u(x, t) is the fluid velocity at point x and time t, 
incomprcssibility is enforced by V.u = 0, Dt = 9* -I- u. V, 
i> is the kinematic viscosity of the fluid, fi the viscos- 
ity parameter for the solute (FENE-P), Tp the polymer 
relaxation time, p the solvent density (set to 1), p the 
pressure, (Vu)^ the transpose of (Vu), Ca/s = (RaRp) 
the elements of the polymer-conformation tensor C (an- 
gular brackets indicate an average over polymer con- 
figurations), T the identity tensor with elements Sap, 
f{rp) = (L^ _ 3)/(L2 - r|,) the FENE-P potential that 
ensures finite extensibility, rp = ^jTr{C) and L the 
length and the maximum possible extension, respectively, 
of the polymers, and c = p,/{v -\- fx) a dimensionless mea- 
sure of the polymer concentration [i^l- c = 0.1 corre- 
sponds, roughly, to lOOppm for polyethylene oxide 0. 

We consider homogeneous, isotropic, turbulence, so we 
use periodic boundary conditions and solve Eg. IT]) by 
using a massively parallel pseudospectral code |23| with 
collocation points in a cubic domain (side L = 27r). 
We eliminate aliasing errors [11 by the 2/3 rule, to 
obtain reliable data at small length scales, and use a 
second-order, slaved Adams-Bashforth scheme for time 
marching. For Eq. we use an explicit sixth-order 
central-finite-difference scheme in space and a second- 
order Adams-Bashforth method for temporal evolution. 
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The numerical error in rp must be controlled by choosing 
a small time step St, otherwise rp can become larger than 
L, which leads to a numerical instability; this time step 
is much smaller than what is necessary for a pseudospec- 
tral DNS of the NS equation alone. Table lists the 
parameters we use. We preserve the symmetric-positive- 
definite (SPD) nature of C at all times by usingj^^ the 
following Cholesky-decomposition scheme: If we define 
J = f{rp)C, Eq. lO becomes 

DtJ = J.(Vu) + {Vuf.J - s{J -I)+ qj, (3) 

where s ^ (L^ + f)/iTpL^). q ^ [d/{L^ - 3) - [L^ - 
S + - 3)/(rpL2(L2 - 3))], f ^ Tr{J), and d - 

Tr[J.(Vu) + iyuY-J]. Since C and hence J are SPD 
matrices, we can write J = CC^ , where £ is a lower- 
triangular matrix with elements £ij , such that £ij = for 



j > i. Thus Eq.® now yields (1 < i < 3 and Ti 
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The SPD nature of C is preserved by Eq.(gj) if > 0, 
which we enforce explicitly |22| by considering the evolu- 
tion of ln(£ii) instead of (a. 

We use the following initial conditions (superscript 
0): C°j„(x) = (5,„„ for aU x; and u°„(k) = 
P„„(k)w°(k)exp(i6'„(k)), with m,n = x,y,z, P„„ = 
{Smn — kmkn/k^) the trausversc projection operator, k 
the wave-vector with components k,n = —N/2 + 

1, . . . ,iV/2), k = |k|, 6'„(k) random numbers distributed 
uniformly between and 27r, and wJJ(k) chosen such 
that the initial kinetic- energy spectra are either of type 
/, with E'{k) = fc2exp(-2fc'*), or of type //, with 
E^'{k) = fc4cxp(-2fc2). 

In addition to u(x, t), its Fourier transform Uk(t), and 
C(x, t) we monitor the vorticity a> = V x u, the kinetic- 
energy spectrum P(fc, i) = Efc-i/2<fc'<fe+i/2 l^k'WI' the 



total kinetic energy £{t) = J2k^(^'^)^ energy- 
dissipation-rate e{t) = v^j^k'^E{k,t), the cumulative 
probability distribution of scaled polymer extensions 
P<^(r|,/L2)^ and the hyperflatness T&{r) = SQ{r)/Sl{r), 
where Sp{r) = ({[u(x -t- r) — u(x)] • r/r}P) is the order-p 
longitudinal velocity structure function and the angular 
brackets denote an average over our simulation domain 
at tm- For notational convenience, we do not display the 
dependence on c explicitly. 

Figure (^) shows that e first increases with time, 
reaches a peak, and then decreases; for c = this peak 
occurs aX t = t,n- The position of this peak changes 
mildly with c but its height goes down significantly as 
c increases. This suggests the following natural defini- 
tion 01 of the percentage drag or dissipation reduction 
for decaying homogeneous, isotropic turbulence: 



DR = 



X 100; 



(5) 



here (and henceforth) the superscripts / and p stand, re- 
spectively, for the fluid without and with polymers and 
the superscript m indicates the time tm- Figure ^p) 
shows plots of DR versus c, for the Weissenberg num- 
ber We = Tp^Jef^^ /v ~ 0.35, and versus We, for 
c=l/ll~0.1. DR increases with c in qualitative accord 
with experiments on channel flows (where DR is defined 
via a normalized pressure difference) ; but it drops gently 
as We increases, in contrast to the behavior seen in chan- 
nel flows (in which Tp is varied by changing the polymer). 

Ill decaying turbulence, the total kinetic energy £{t) 
of the fluid falls as t increases; the rate at which it falls 
increases with c [Fig. QJ:)], which suggests that the ad- 
dition of polymers increases the effective viscosity of the 
solution. This is not at odds with the decrease of e with 
increasing c since the effective viscosity because of poly- 
mers turns out to be scale- dependent. We confirm this 
by obtaining the kinetic-energy spectrum i?^''™(fc) for the 
fluid in the presence of polymers at i = tm- For small 
concentrations (c ~ 0.1) the spectra with and without 
polymers differ substantially only in the deep dissipa- 
tion range, where Ef^'^{k) < £'P'™(fc). As c increases, 
to say c ~ 0.4, E^'^^ik) is reduced relative to Ef'"^{k) 
at intermediate values of k [Fig. 13^)]; however, deep in 
the dissipation range Ef^"^{k) < E^'^^ik). We now de- 
fine the effective scale-dependent viscosity Ve{k) = 

V + Avlk), with Al/(fc) = -tJ'J2k-l/2<k'<k+l/2 Uk' ■ (V ■ 

J)-]^'/[Tpk'^EP'"'{k')], where (V • J7)k is the Fourier 
transform of V • J". The inset of Fig. shows that 
Aj/(fc) > for fc < 15, but Ai^(fc) < around k = 
20. This explains why EP''^{k) is suppressed relative to 
Ef'^^ik) at small fc, rises above it in the deep-dissipation 
range, and crosses over from its small-fc to large-fc be- 
haviors around the value of fc where Ai'(fc) goes through 
zero. 

Given the resolution of our DNS, inertial-range inter- 
mittency can be studied only by using extended self sim- 
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ilarity [2^ as we will report elsewhere. However, we 
explore dissipation-range statistics further by calculat- 
ing the hyperflatness ^^{r) [Fig. (jSt)]- The addition of 
polymers slows down the growth of J-Q{r)^ as r —* 0, 
which signals the reduction of small-scale intermittency. 
This is further supported by the iso-|a;| surfaces shown 
in Fig. (PJ. If no polymers are present, these iso-|w| sur- 
faces are filamentary [i^ for large \uj\; polymers suppress 
a significant fraction of these filaments. 

We use a rank-order method [2^ to obtain P'-^(rp/L^) 
and find that, as c increases [Fig. (0:)], the extension of 
the polymers decreases. We have checked that, in the 
passive-polymer version of Eqs.(^ and the extension 
of polymers is much more than in Fig. ||21:). 

Our study contrasts clearly drag reduction in homoge- 
neous, isotropic, turbulence and in wall-bounded flows. 
In both these cases the polymers increase the overall vis- 
cosity of the solution (see, e.g., Fig. 1^;) and Ref.0|). 
In wall-bounded flows the presence of polymers inhibits 
the flow of the stream-wise component of the momentum 
into the wall, which, in turn, increases the net through- 
put of the fluid and thus results in drag reduction, a 
mechanism that can have no analog in homogeneous, 
isotropic turbulence. However, the decrease of e{t) with 
increasing c [Fig. lH))] yields a natural definition of DR 
[Eq.|5l)] for this case [27|. Thus, if the term drag reduc- 
tion must be reserved for wall-bounded flows, then we 
suggest the expression dissipation reduction for homo- 
geneous, isotropic, turbulence. We have shown that 
must be scale-dependent; its counterpart in wall-bounded 
flows is the position-dependent viscosity of Rcfs. 0,0- 
Furthermore, as in wall-bounded flows, an increase in c 
leads to an increase in DR [Fig. (^)]. In channel flows 
an increase in We leads to an increase in DR, but we find 
that DR falls marginally as We increases [Fig. ^p)]. 

Our DNS of the Navier-Stokes equation with poly- 
mer additives [Eqs. and (jSjl] resolves the controversy 
about drag reduction in decaying homogeneous, isotropic 
turbulence and shows clearly that Eq. Q offers a natu- 
ral definition of DR for this case in a far more realistic 
model than those of Refs. [ill H^. We also find a non- 
trivial modification of the fluid kinetic-energy spectrum 
especially in the deep-dissipation range [Fig. ^fl)] that 
can be explained in terms of a polymer- induced, scale- 
dependent viscosity. Experiments |l6l llTj do not resolve 
the dissipation range as clearly as we do, so the experi- 
mental verification of the deep-dissipation-range behav- 
ior of Fig. i B b, ) remains a challenge. Earlier theoretical 
studies [iJJ, |ll| have also not concentrated on this dis- 
sipation range. The reduction in the small-scale inter- 
mittency [Fig. (03)] and in the constant- |a;| isosurfaces 
[Fig. (Pll is in qualitative agreement with channel-flow 
studies j3| , where a decrease in the turbulent volume frac- 
tion is seen on the addition of the polymers, and water-jet 
studies [iJl, where the addition of the polymers leads to 
a decrease in small-scale structures. We hope our work 
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TABLE I: The parameters A^, St, L, v, rp and c for our 
four runs NSP - 96, NSP - 192, NSP - 256A, and NSP - 256. 
NSP - 96, NSP - 192, NSP - 256A use type / initial conditions; 
NSP — 256 uses an initial condition of type //. We also carry 
out DNS studies of the NS equation with the same nu merical 
resolutions as our NSP runs. 

and We = Tp^ef^^jv; NSP-96: Re = 47.1 and We = 
0.03, 0.17, 0.24, 0.28, 0.31, 0.41, 0.48, 0.55, 0.62, 0.68, 1.03; 
NSP-192 and NSP-256A: Re = 47.1 and We = 0.35; NSP-256: 
Re = 126.6 and We = 0.76. 



will stimulate more experimental studies of drag or dis- 
sipation reduction in homogeneous, isotropic turbulence. 
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